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Abstract. We demonstrate the existence of stochastic resonance (SR) in confined systems arising from 
entropy variations associated to the presence of irregular boundaries. When the motion of a Brownian 
particle is constrained to a region with uneven boundaries, the presence of a periodic input may give rise 
to a peak in the spectral amplification factor and therefore to the appearance of the SR phenomenon. We 
have proved that the amplification factor depends on the shape of the region through which the particle 
moves and that by adjusting its characteristic geometric parameters one may optimize the response of 
the system. The situation in which the appearance of such entropic stochastic resonance (ESR) occurs 
is common for small-scale systems in which confinement and noise play an prominent role. The novel 
mechanism found could thus constitute an important tool for the characterization of these systems and 
can put to use for controlling their basic properties. 

PACS. 02.50.Ey Stochastic processes - 05. 40. -a Brownian motion - 05.10. Gg Stochastic analysis methods 



1 Introduction 



A Brownian particle moving in a threshold-like potential 
landscape and subjected to the influence of a periodic forc- 
ing may exhibit a coherent response giving rise to an am- 
plification of the input at a certain optimal value of the 
noise level. This resonant phenomenon, observed in gen- 
eral in the wide class of periodically modulated noisy sys- 
tems, was termed stochastic resonance and constituted a 
paradigm shift in the way we think about noise effects in 
systems away from equilibrium [1 . In this new paradigm, 
the presence of noise does not always constitute a nui- 
sance; on the contrary, it may play a constructive role 



Up to now, the phenomenon of SR has been observed 
mainly in systems dominated by the presence of a purely 
energetic potential or possessing some dynamical thresh- 
old pQ. However, when scaling down the size of a sys- 
tem, the free energy rather than the internal energy be- 
comes the most appropriate potential, and there are cases 
in which changes in the free energy are mainly due to 
entropy variations [11113114115116117118] . This is what oc- 
curs in constrained systems. In the case of a Brownian 
particle moving in a confined medium, entropy variations 
contribute to changes in the free energy and may un- 
der some circumstances become its leading contribution 
|12|13|17|18] , We will show in this work that the uneven- 
ness may also give rise to a stochastic resonance effect 
and that this effect can be controlled upon variation of 



the geometrical parameters, characterizing the shape of 
the cavity in which the Brownian particle dwells. 

Usually, the analysis of SR effects have been performed 
by means of pertinent Langevin or corresponding Fokker- 
Planck models [19120) . In confined systems the presence of 
boundaries exerts a strong influence in the dynamics and 
one has to solve the corresponding boundary value prob- 
lem. This task cannot always be easily achieved. The fact 
that in many instances boundaries are very intricate enor- 
mously complicates the mathematical treatment of the 
problem to the extent of becoming a Herculean task when 
the boundaries are extremely irregular [21] . This feature 
demands the implementation of different approaches en- 
tailing a simplification of the analysis 22 23 24J . Among 
them, the Fick- Jacobs equation, based on a coarsening of 
the description in terms of a single, relevant coordinate de- 
gree of freedom, accurately performs this task [17|18|21|22|23j . 
This methodology will guide us in this article to analyze 
the appearance of the SR effect in presence of unevenness. 

The article is organized in the following way. In Sec- 
tion 2, we introduce a model for Brownian motion in a 
confined medium. In Section 3, we present a reduction 
method which simplifies the complex nature of the 3D/2D 
dynamics giving rise to an effective one-dimensional ki- 
netic description. Section 4 is devoted to evaluate the tran- 
sition rate from the reduced kinetic description, and the 
introduction of spectral amplification within the two-state 
approximation. In Section 5, we present the results show- 
ing the ESR phenomenon. The impact of the geometrical 
shape and confinement on the spectral amplification is 
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Fig. 1. Schematic illustration of the two-dimensional structure 
confining the motion of the Brownian particles. The symmet- 
ric structure is defined by a quartic double well function, cf. 
Eq. ([2]), involving the geometrical parameters L x , L y and b. 
Brownian particles are driven by a sinusoidal force F(t) along 
the longitudinal direction and a constant force G in the or- 
thogonal direction. 

discussed in Section 6. Finally, we summarize our main 
conclusions in Section 7. 



2 Confined Brownian motion 



very close to the lower boundary of the structure, recov- 
ering the effect of an energetic bistable potential. For the 
sake of a dimensionless description, we measure all lengths 
in units of L x , i.e. x — x/L x , y = y/L x implying b = b/L x 
and w\ = w\/L x = —w u , temperature in units of an ar- 
bitrary, but irrelevant reference temperature Tr and time 
in units of r = yL x /k^T^ that is, twice the time the par- 
ticle takes to diffuse a distance L x at temperature Tr, i.e. 
i = t/r and fl = fir. We scale forces by Fr = jL x /t, i.e. 
the orthogonal force reads G = G/Fr and the sinusoidal 
force F(i) = F(t)/F-R. For better legibility, we shall omit 
the tilde symbols in the following. In dimensionless form 
the Langevin-equation {1} and the boundary functions ([2]) 
read: 

dr 

— = -Ge y - F(t)e x + , (3) 

w\(x) = -w u (x) = ex 4 - 2ex 2 - b/2 , (4) 

where we defined the aspect ratio e = L y /L x and the 
rescaled temperature D = T/Tr. 



The dynamics of a particle in a constrained geometry sub- 
jected to a sinusoidal oscillating force F(t) along the axis 
of the structure and to a constant force G acting along 
the orthogonal, or transverse, direction can be described 
by means of the Langevin equation written, in the over- 
damped limit, as 



7 



dr 

dt 



Ge y - F(t)e x + y/ykaT^t) , 



(1) 



where r denotes the position of the particle, 7 is the fric- 
tion coefficient, e x and e y the unit vectors along x and y- 
directions, respectively, and £(t) is a Gaussian white noise 
with zero mean which obeys the fluctuation-dissipation 
relation (&(*)£,•(*')) = 2 for ij = x,y. The 
explicit form of the longitudinal force is given by F(t) = 
F sin(nt) where F is the amplitude and Q is the angular 
frequency of the sinusoidal driving. 

In the presence of confinement, this equation has to 
be solved by imposing reflecting (no-flow) boundary con- 
ditions at the walls of the structure. For the 2D structure 
depicted in Fig. [TJ the walls are defined by 



W\{x) = Ly 



L, 



2L„ 



-w u (x), (2) 



where w\ and w u correspond to the lower and upper bound- 
ary functions, respectively. The characteristic length L x 
refers to the distance between the bottleneck and the po- 
sition of maximal width, L y corresponds to the narrowing 
of the boundary functions and b to the remaining width at 
the bottleneck, cf. Fig.[TJ Consequently, the local width of 
the structure reads: 2w(x) = w u (x) — w\(x) . This particu- 
lar choice of the geometry is intended to resemble the clas- 
sical setup for SR in the context of energetic barriers. In 
fact, in the limit of a sufficiently large transverse force G, 
the particle is in practice restricted to explore the region 



3 Reduction of dimensionality 

Since the above mentioned dynamics given by Eq. ([3|) with 
the boundary conditions could not be solved analytically, 
we simplified the problem by assuming equilibration in y- 
direction and thereby reducing the dimensionality of the 
problem |22l23l24l25j . 

First, we consider the case in the absence of the peri- 
odic forcing, i.e. F(t) = 0. Then, the 2D dynamics is de- 
scribed by the following 2D Smoluchowski equation [20|19j 

|p(x, y, t) = ^e-^/^e^^/^Pfe y, t) 

+ D^-e- u ^l D ^e u ^y D P{x, y, t) , (5) 

with reflecting boundary conditions at the confining walls 
and where the potential function is given by U{x,y) = 
G y. Since we are mainly interested in the dynamics in re- 
direction, we introduce the marginal probability density 
P{x, t) which is obtained by integration over the trans- 
verse coordinate: 



(6) 



P(x,t) = / dyP(x,y,t). 



On integrating Eq. ([5]) over the transverse direction, we 
get 

|P(M) = 

D lL J dy {z- u[x ' v),D ^ u{x > y)/D P^y,t)} . (7) 

Assuming local equilibrium in the y-direction, we define 
the ^-dependent effective energy function A(x) (omitting 
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irrelevant constants) reading 

-A(x)/D = f dye -U(x, y )/D 



(8) 



Consequently, the conditional local equilibrium probabil- 
ity distribution of y at a given x becomes 



p{y; x) 



e -U(x,y)/D e A(x)/D 



(9) 



and is normalized for every i. As a result, the two di- 
mensional probability distribution can be approximately 
expressed as 



P(x,y,t) ^ P(x,t)p(y;x) , 



(10) 



and the kinetic equation for the marginal probability dis- 
tribution, cf. Eq. |(7|), becomes 

I^o-d^'^'^o. (ID 

In the present case with a constant force G in the negative 
y-direction the potential function A(x) reads, cf. Eq. ([8]): 



A(x) 



Din 



-Gy/D 



dy 



= -D In 



wi(x) 

D 
G 



e -Gwi{x)/D _ e -Gw u (x)/D 



(12) 



Making use of the symmetry of our considered structure, 
i.e. w u (x) = —w\(x) and the definition of the half width 
function w(x) = (u> u — W\)/2, the potential function turns 
into 



A(x) 



-D In 



2D . 
— smh 



Gw(x) 
D 



Eq. pTjl can then be rewritten as 
dP(x,t) 



dt 



(13) 



(14) 



with the potential function A(x) given by Eq. (fT2"|) (or by 
Eq. (| 1 3 p for w\(x) = —w u (x)) and with the prime refer- 
ring to the derivative with respect to x. In general, after 
the coarse-graining the diffusion coefficient will depend on 
the coordinate x, but since in our case <C 1, the 

correction can be safely neglected |21l23l24l25l26j . 

It is important to highlight that the potential A(x) was 
not present in the 2D Langevin dynamics, but arises due 
to the entropic restrictions associated to the confinement. 
Then, equation (fT4|) describes the motion of a Brownian 
particle in a free energy potential of entropic nature, as 
A(x) does not only depend on the energetic contribution of 
the force G, but also on the temperature D and the geom- 
etry of the structure in a non-trivial way. For a structure 
like that depicted in Fig. [1] the free energy A(x) forms a 
double-well potential, cf. Fig. [51 As the width at the bot- 
tleneck of the channel approaches zero, the potential A(x) 
diverges. 

It is worth to analyze the two limiting situations that 
can be obtained depending on the transverse force G. 
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Fig. 2. (Color online) The effective ID potential A(x) is de- 
picted (solid black line) for G = 1 and for different noise 
strengths. For comparison the potential functions for the pure 
energetic case (dotted blue line), i.e. A(x) — Gw\(x), and 
for the pure entropic case (dashed red line), i.e. A(x) = 
— Dln[2w(x)], are plotted as well. 



Energy-dominated situation: For G 3> 1, Eq. (fT3|) 

yields A(x) — Gw\(x) (neglecting irrelevant constants) 
and the ID kinetic equation (| 14[) becomes the stan- 
dard Fokker-Planck equation for Brownian motion in 
a purely energetic potential whose shape resembles the 
lower boundary of the structure, 



dP(x,t) 
dt 



d ( dP(x,t) 
dx | dx 



Gw[(x)P(x,t) } . (15) 



Entropy-dominated situation: In the opposite limit, 
i.e. for G <C 1, the effective potential is dominated by 
the purely entropic contribution A{x) = —D ln[2 w(cc)] 
and the kinetic equation turns into the Fick- Jacobs 
equation [22], 



dP(x,t) 
dt 



d_ 

dx 



D 



dP{x,t) 
dx 



(x) „, 
D—^P(x,t) 



(16) 



4 Two-State approximation 

It is instructive to analyze the occurrence of stochastic res- 
onance in the context of the two-state approximation |27j . 
Accordingly, the ID kinetics given by Eq. (fTf)) can be ap- 
proximately mapped into a two-state system with the two 
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Fig. 3. (Color online) The transition rate for the escape from 
one basin to the other obtained after applying a Two-State 
approximation of the ID problem, cf. Eq. (| lTf) and Eq. (|21l) . 
is depicted as a function of the thermal noise D for G = 1.0. 
The 2D structure is denned by the boundary function given 
in Eq. (|4j with e = 1/4 and b = 0.02. Here, the Kramers- 
Smoluchowski rate (K), i.e. tk, converges to the rate given by 
the inverse mean first passage time (MFPT) in the low noise 
limit. 



states corresponding to the two wells of the symmetric ef- 
fective potential A(x). An estimate for the transition rates 
could be obtained by applying the Mean-First-Passage- 
Time (MFPT) approach [25] • 



4.1 Mean First Passage Time approach 



In order to calculate the transition rate from one state to 
the other, one evaluates the inverse of the MFPT to reach 
a potential minimum after starting out from the other 
minimum of the symmetric and bistable potential A(x). 
Then, the transition rate is given by 

r MF P T (A G) = Tjrr^—f) ' (17) 

where T\(— 1 — * 1) is the first moment of the first passage 
time distribution for reaching x = 1 starting out at x = 
— 1. The nth moment of the first passage time distribution 
obeys the following recurrence relation [28] 

T n (-1 -> 1) := (t"(-l -> 1)) 

= " f Axe A ^' D [ X dye- A ^/° 

V J-l J-oo 

x(t n -\y^l)) (18) 

for n € N and with T (a — > b) — 1 for arbitrary a and b. 
Accordingly, within the one dimensional approximation, 
cf. Eq. (fl4|) . the mean first passage time (i.e., n = 1) for 



the potential function A(x), cf Eq. (|T3|) . reads 

(19) 

where x\ is the left limiting value at which the bound- 
ary function vanishes. Eq. (fl9|) can be evaluated using a 
steepest descent approximation leading to the commonly 
known Kramers- Smoluchowski rate. 



4.2 Kramers-Smoluchowski rate 

For a potential A(x) with a barrier height A A ^> D 
the escape rate of an overdamped Brownian particle from 
one well to the other in the presence of thermal noise, 
and in the absence of a force, is given by the Kramers- 
Smoluchowski rate [27 28 29 30 . reading in dimcnsionless 
units, 

rm VA"(x min )\A"(x m ^)\ f —AA\ 
MD) = ^ exp ^— j , (20) 

where A" is the second derivative of the effective potential 
function with respect to x, and x max and x m i n indicate the 
position of the maximum and minimum of the symmetric 
potential, respectively. 

For the potential given by Eq. (fT3|) and the shape de- 
fined by Eq. ([¥]), the corresponding Kramers-Smoluchowski 
rate for transitions from one basin to the other reads [181 



r K (D) = 




Note that the Kramers-Smoluchowski approximation 
yields good results for barrier heights AA much larger 
than the thermal energy that in the present scaling is given 
by D. In fact, for A A ^> D the Kramers-Smoluchowski 
rate rx approximates accurately the rate tmfpt evaluated 
numerically from the MFPT expression ([T9|) . as depicted 
in Fig. H 



5 Role of the transverse force G 

In section [3] we introduced the Fick-Jacobs equation to 
approximatively describe the Brownian motion in a 2D 
structure like that depicted in Fig. Q] using a simplified ID 
modeling with an effective bistable potential. This poten- 
tial exhibits a barrier the particle has to overcome noisily 
in order to make a transition from one well to the other. 
For a sinusoidal driving force applied in .T-direction, i.e. 
F(t) — Fq sm(f2t), a synchronization effect between the 
oscillatory forcing and the noise-induced transitions over 
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directly from the numerical integration of the ID kinetic 
equation. In the presence of an oscillating force F(t) in 
a;— direction there is an additional contribution to the ef- 
fective potential in Eq. (fT4|) and the ID kinetic equation 
in dimensionless units reads 1181 



Fig. 4. (Color online) A detailed comparison of the behav- 
ior of the spectral amplification factor obtained numerically 
from the ID modeling (solid line), cf. and Eq. (f25)) , 

and within the two-state approximation (Eq. (|22[) b using either 
the Kramers-Smoluchowski rate (K) (dotted line) or the rate 
obtained directly from the mean first passage time (MFPT) 
(dashed-line) . Here, the transverse force is G — 1.0, and the 
structure is defined by Eq. Q with an aspect ratio e = 1/4, 
and a bottleneck width b = 0.02. 



the entropic potential barrier takes place and was reported 
previously in Ref. [18 . Under these circumstances and in 
the presence of a finite orthogonal, i.e. transverse force G, 
increasing the noise level results in a noise-amplified re- 
sponse signal. The improvement of the response is quanti- 
fied in the following by the spectral amplification factor r), 
which is the ratio of the power stored in the response of the 
system at the driving frequency fl and the power of the 
sinusoidal driving signal. The occurrence of the Entropic 
Stochastic Resonance effect manifests in the presence of an 
optimal dose of noise for which the spectral amplification 
is maximal 18 . 



5.1 Two-State modeling 

It is straight forward to derive an analytic expression for 
the spectral amplification within a two-state modeling. 
The sinusoidal driving modulates the transition rates, which 
are given either by tmfpt or tk (cf . Sec. [4J . Within a first 
order perturbation theory in the ratio of driving ampli- 
tude Fq and noise level D, it is possible to find a closed 
expression for the response of the two-state system and ac- 
cordingly for the spectral amplification factor that reads 



1 



v 



4r| ate (£>) 



D 2 4^(25) 



n 2 ' 



(22) 



where rR ato is the transition rate of the unperturbed two- 
state system and is given by Eq. (fTF]) , or approximately 
by Eq. (1201. 



5.2 ID modeling 

Avoiding the approximations involved in the two-state 
modelling, the system's response could also be obtained 



(23) 



By spatial discretization, using a Chebyshev colloca- 
tion method, and employing the method of lines, we re- 
duced the kinetic equation to a system of ordinary differ- 
ential equations, which was then solved using a backward 
differentiation formula method [31] • As a result, we ob- 
tained the time dependent probability distribution P(x, t) 
and, from that, the time-dependent average position de- 
fined as 



(*(*)> 



x P{x, t)dx . 



(24) 



which was computed in the long-time limit. After a Fourier- 
expansion of (x(t)) one finds the amplitude Mi of the 
first harmonic of the system's response. Hence, the spec- 
tral amplification r\ [30] for the fundamental oscillation is 
evaluated as 



V 



Mi 



(25) 



The spectral amplification r\ depicts a bell-shaped be- 
havior, cf. Fig. [4] indicating the existence of a Stochastic 
Resonance effect: there is a maximum of the spectral am- 
plification at an optimal value of noise. The qualitative 
behavior is captured by the two-state approximation as 
well. The prediction of the two-state modelling with rates 
evaluated directly from MFPT is closer to the ID mod- 
elling than those obtained by making use of the Kramers- 
Smoluchowski approximated rate, as illustrated in Fig. 0] 
This is due to the fact that the condition AA 3> D is 
not always fulfilled in this case. In fact for the entropy- 
dominated situation, AA and D are of the same order of 
magnitude. 

Next we present the results of numerical simulations of 
the full (2D) problem. By doing so we demonstrate that 
the ESR effect is robust and not just an artifact of the 
reduction procedure. 



5.3 2D modelling 

The accuracy of the reduced one-dimensional kinetic de- 
scription can be examined by comparing the results with 
those obtained by Brownian dynamic simulations, per- 
formed by integration of the overdamped 2D Langevin 
equation (|TJ) . The simulations were carried out using the 
standard stochastic Euler-algorithm. 

The resulting amplification factor as a function of the 
value of the transverse force and the noise strength are 
plotted in Fig. [5] For a finite transverse force G the spec- 
tral amplification exhibits a peak at an optimum value of 
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Fig. 5. (Color online) The spectral amplification obtained 
from 2D simulations is depicted as a function of both the trans- 
verse force G and the noise strength D for an input signal fre- 
quency Q — 0.1 and amplitude Fq = 0.1. The shape of the 
two-dimensional channel is defined by the dimensionless func- 
tion w(x) = —ex 4 + 2 ex' 2 + b/2 with the aspect ratio e = 1/4 
and bottleneck width b = 0.02. The symbols connected by 
lines correspond to the same G-values. In particular, these are 
G = 0.0, 0.05, 0.1, 0.125, 0.15 and 0.2. In the deterministic 
limit, i.e D — » 0, the spectral amplification reaches a limit- 
ing value (I/O 2 ) for G = 0. For G = the effect of Entropic 
Stochastic Resonance disappears. 



the noise strength which is indicative of the effect of En- 
tropic Stochastic Resonance. However, for vanishing trans- 
verse force G the spectral amplification does not exhibit 
any peak and decays monotonically with increasing noise 
level, i.e. the ESR - effect is not observed. In the determin- 
istic limit, i.e. D — ► 0, the spectral amplification reaches 
the limit value r\ — 1/ Q 2 for G = 0, as shown in Fig. [5l 
We remark that the later is only true when the amplitude 
of the system's response, which is the ratio of input signal 
amplitude to input signal frequency, is smaller than the 
value xi, which is the limiting value at which the bound- 
ary function vanishes. It is also worth to point out that 
starting out from a finite G-value the position of the ESR 
peak shifts towards smaller noise strengths as G decreases 
while the maximum value increases, cf. Fig. A compar- 
ison between the ID-modelling and the full 2D simulation 
is depicted in Fig. [6] The 2D simulation results convinc- 
ingly corroborate the validity of the modelling within the 
ID Fick-Jacob approximation, see also below. 

Thus, we detect a non-monotonic behavior of the spec- 
tral amplification only for finite values of the orthogonal 
force G/0, while for G = the spectral amplification de- 
cays monotonically. In other words, for the dynamics and 
situation considered in Sec. [2] the occurrence of the ESR- 
effect requires a non- vanishing orthogonal force G which 
contributes within the ID modelling to the effective en- 
tropic potential. 



D 

Fig. 6. (color online) The dependence of the spectral amplifi- 
cation n on noise level D, at two different values of the trans- 
verse force G, at a constant input frequency and amplitude. 
The symbols correspond to the results of the Langevin simu- 
lations for the two-dimensional channel whose shape is defined 
by w(x) = —ex 4 + 2 ex 2 + 0.01 with an aspect ratio e — 1/4, 
whereas the lines are the results of the numerical integration 
of the ID kinetic equation (|23jl . In the deterministic limit, i.e 
D — » 0, the spectral amplification reaches the limiting value 
(l/« 2 ) for G = 0. 




x 



Fig. 7. Schematic illustration of the two-dimensional structure 
confining the motion of the Brownian particles. The shape de- 
fined by the dimensionless function w(x) = 12. 5x s — 27.5a; 6 + 
18.21s 4 - 3.92x 2 - 0.01. 

6 Role of the shape of the structure 

Fig. [6] depicts the main findings for the occurrence of the 
ESR-effect in an effective bistable configuration with two 
large basins connected by a small bottleneck, cf. Fig. [1] 
Namely: (i) the occurrence of the ESR-effect due to the 
interplay of a finite value of the orthogonal force G, and 
(ii) the ability of the ID Fick- Jacobs approximation to re- 
produce very accurately the results of the full 2D problem 

mm- 

Next, in order to investigate the impact of the shape of 
the structure on the ESR behavior we have considered yet 
another channel geometry: It is similar to the one depicted 
in Fig. [1] but exhibits two intermediate wider bottlenecks 
which are connected via a much narrower bottleneck, see 
in Fig. The geometric shape of the structure is defined 
by the dimensionless width function 

w(x) = 12.5a; 8 - 27.5a; 6 + 18.21a; 4 - 3.92x 2 - 0.01 . (26) 

Following the same analysis performed in Sec.O we ob- 
tained the spectral amplification within the ID modelling 
and compared it with results obtained from 2D numerical 
simulation of the Langevin equation |3]) in the new 2D 
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Fig. 8. (Color online) The dependence of the spectral am- 
plification r\ on noise level D, at two different values of the 
transverse force G, at a constant input frequency and am- 
plitude for the channel illustrated in Fig. [7J. The symbols 
correspond to the results of the Langevin simulations for the 
two-dimensional structure whose shape is defined by w(x) — 
12. 5x 8 - 27.5s 6 + 18.21x 4 - 3.92s 2 - 0.01, whereas the lines are 
the results of the numerical integration of the ID kinetic equa- 
tion (|23[) . Like in Fig. [6] in the deterministic limit (D — > 0), 
the spectral amplification reaches the limiting value (1/ ' Q 2 ) for 
G = 0. 



structure defined by Eq. (f2"o| . Again, the effective poten- 
tial function A(x) obtained from the Fick- Jacobs approx- 
imation exhibits a large potential barrier separating the 
two basins each of which is additionally separated by a 
smaller potential barrier into two wells. In fact, the con- 
struction of the alternative geometric structure was done 
in such a way to reveal the existence of a main entropic 
barrier controlling the transitions between the two main 
basins. 

The behavior of spectral amplification as a function of 
the noise strength for the new channel defined by Eq. 1261 
and for two different values of the transverse force is de- 
picted in Fig. [8] Interestingly, we can observe the existence 
of Entropic Stochastic Resonance even in the new struc- 
ture, when the transverse force is present in the system. 
Comparing the results of Fig. [8] for the channel depicted 
in Fig. [7] with those in Fig. [|)J one can observe that the 
ESR peak appears at higher values of the noise strength. 
In addition, the enhancement of the signal upon increas- 
ing the noise is more pronounced in the new structure. 
These results suggest that by a proper design of the ge- 
ometry of the channel it would be possible to significantly 
enhance and optimize the response of a confined system. 
This would be specially important in biological systems, 
where the noise (i.e. the temperature) is a variable that 
can neither be arbitrarily chosen nor eliminated. Finally, 
it is also worth stressing that the Fick-Jacobs approxima- 
tion still holds nicely in this case since the 2D numerical 
simulation results are in very good agreement with those 
obtained by ID modelling. 

Overall, the results indicate that the existence of a 
small bottleneck separating two basins and forming an 
effective entropic potential barrier, leads robustly to the 
occurrence of an Entropic Stochastic Resonance effect. 



7 Conclusions 

We have shown that unevenness may be the origin of 
many resonant phenomena in small-scale systems. The 
constrained motion of a Brownian particle in a region 
limited by irregular boundaries impedes the access of the 
particle to certain regions of space giving rise to entropic 
effects that can effectively control the dynamics. The in- 
terplay between the noise present in the system, the ex- 
ternal modulation and the entropic effects results in an 
entropic stochastic resonance. However, the presence of a 
transverse force G is crucial to observe this resonant be- 
havior. This ESR is genuine of small scale systems where 
confinement yields entropic effects. The occurrence of ESR 
depends on the shape of the channel and can then be con- 
trolled by it. Thus understanding the role of noise and 
confinement in these systems does provide the possibil- 
ity for a design of stylized channels wherein response and 
transport become efficiently optimized. 

This work has been supported by the DFG via research 
center, SFB-486, project A10, the Volkswagen Founda- 
tion (project 1/80424), the German Excellence Initiative 
via the Nanosy 'stems Initiative Munich (NIM), and the 
DGCyT of the Spanish government through grant No. 
FIS2005-01299. 



References 

1. L. Gammaitoni, P. Hanggi, P. Jung, F. Marchesoni, Rev. 
Mod. Phys. 70, 223 (1998) 

2. A. Bulsara, P. Hanggi, F. Marchesoni, F. Moss, M. 
Shlesinger, J. Stat. Phys. 70, 1 (1993) 

3. A.R. Bulsara, L. Gammaitoni, Phys. Today 49 (3), 39 
(1996) 

4. P. Hanggi, ChemPhysChem 3, 285 (2002) 

5. J.M.G. Vilar, J.M. Rubf, Phys. Rev. Lett. 77, 2863 (1996); 
J.M.G. Vilar, J.M. Rubf, Phys. Rev. Lett. 78, 2886 (1997) 

6. V.S. Anishchenko, A.B. Neiman, F. Moss, L. Schimansky- 
Geier, Phys. Usp. 42, 7 (1999) 

7. G. Schmid, I. Goychuk, P. Hanggi, Europhys. Lett. 56, 22 
(2001) 

8. T. Wellens, V. Shatokhin, A. Buchleitner, Rep. Prog. Phys. 
67, 45 (2004) 

9. H. Yasuda, T. Miyaoka, J. Horiguchi, A. Yasuda, P. Hanggi, 
Y. Yamamoto, Phys. Rev. Lett. 100, 118103 (2008) 

10. C.R. Nicolis, G. Nicolis, Scholarpedia, 2(11), 1474 (2007) 

11. I. Goychuk, P. Hanggi, Phys. Rev. Lett. 91, 070601 (2003) 

12. I. Goychuk, P. Hanggi, J. L. Vega, S. Miret Artes, Phys. 
Rev. E 71, 061906 (2005) 

13. B. Hille, Ion Channels of Excitable Membranes (Sinauer, 
Sunderland, 2001) 

14. R.M. Barrer, Zeolites an Clay Minerals as Sorbents and 
Molecular Sieves (Academic Press, London, 1978) 

15. L. Liu, P. Li, S.A. Asher, Nature 397, 141 (1999) 

16. A.M. Berezhkovskii, S.M. Bezrukov, Biophys. J. 88, 
L17(2005) 

17. D. Reguera, G. Schmid, P.S. Burada, J.M. Rubf, P. 
Reimann, P. Hanggi, Phys. Rev. Lett. 96, 130603 (2006) 

18. P.S. Burada, G. Schmid, D. Reguera, M.H. Vainstein, J.M. 
Rubi, P. Hanggi, Phys. Rev. Lett. 101, 130602 (2008) 



8 



Burada et al: Entropic Stochastic Resonance 



19. P. Hanggi, H. Thomas, Phys. Rep. 88, 207 (1982) 

20. H. Riskcn, The Fokker-Planck equation, 2nd ed. (Springer, 
Berlin, 1989) 

21. P.S. Burada, G. Schmid, D. Reguera, J.M. Rubf, P. 
Hanggi, Phys. Rev. E 75, 051111 (2007); P.S. Burada, G. 
Schmid, P. Talkner, P. Hanggi, D. Reguera, J.M. Rubf, 
BioSystems 93, 16 (2008) 

22. M. H. Jacobs, Diffusion Processes (Springer, New York, 
1967) 

23. R. Zwanzig, J. Phys. Chem. 96, 3926 (1992) 

24. D. Reguera, J.M. Rubf, Phys. Rev. E 64, 061106 (2001) 

25. P. Kalinay, J.K. Percus, Phys. Rev. E 74, 041203 (2006) 

26. A. M. Berezhkovskii, M.A. Pustovoit, S.M. Bezrukov, J. 
Chem. Phys. 126, 134706 (2007) 

27. B. McNamara, K. Wiesenfeld, Phys. Rev. A 39, 4854 
(1989) 

28. P. Hanggi, P. Talkner, M. Borkovec, Rev. Mod. Phys. 62, 
251 (1990) 

29. H. Kramers, Physica (Utrecht) 7, 284 (1940) 

30. P. Jung, P. Hanggi, Phys. Rev. A 44, 8032 (1991); ibid, 
Europhys. Lett. 8, 505 (1989) 

31. NAG Fortran Library Manual, Mark 20 (The Numerical 
Algorithm Group Limited, Oxford, England, 2001) 



